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Abstract 

We  analyze  the  "zero  crossings"  of  a  solution  to  the  heat 
equation.  Zero  crossings  of  differences  of  multiple  levels  of 
resolution  of  data  have  been  used  for  representation  and  analysis  of 
digitized  signals  and  images.  We  establish  a  mathematical  structure, 
based  on  the  heat  equation,  for  investigating  these  kinds  of 
representations.  We  show  that  zero  crossings  are  surfaces  which  evolve 
in  time,  and  that  no  new  zero  crossing  surfaces  appear  at  positive 
values  of  time.  We  show  that  complete  knowledge  of  the  zero  crossings 
of  a  solution  to  the  heat  equation,  together  with  gradient  data  at 
those  zero  crossings,  uniquely  determines  the  initial  data.  We  provide 
an  algorithm,  via  the  backwards  heat  equation,  for  reconstructing  the 
initial  data  given  the  zero  crossings  and  gradient.  We  argue  that,  in 
general,  the  zero  crossings  alone  are  insufficient  to  determine  the 
initial  data. 
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I.   The  Heat  Equation  for  Multlcesolutlon  RepFeaentatton 

It  is  increasingly  clear  that  the  use  of  multi-resolution 
representation  is  an  important  idea  for  the  analysis  of  signal  and 
image  data.  Many  data  structures  have  been  studied,  including 
"gaussian  pyramids,"  difference  of  gaussian  channels,  and  Laplacian 
pyramids  [1,2,3].  The  goal  of  these  decompositions  of  scalar  data  is 
the  representation  of  essential  information  In  a  scale  invariant 
fashion  permitting  easy  feature  extraction  and  matching  with  models. 
We  show  below  that  these  structures  may  be  studied  by  means  of  a 
unifying  mathematical  treatment  based  on  the  heat  equation.  We  then 
define  the  notion  of  zero-crossings,  and  discuss  their  significance. 

Let  f(x)  be  a  bounded  function  defined  for  xSlR".  If  we  blur 
f(x)  by  increasingly  diffuse  gaussians  parameterized  by  t  >  0,  we 
obtain 

u(x,t)  =  j    K(x-y,t)f(y)dy  ,  (1.1) 

IR" 

where 

K(x.t)  =  ('iTrt)-n/2  e-lx|2/4t  .  (1.2) 

K(x,t)  is  the  fundamental  solution  to  the  heat  equation;  in  fact, 
u(x,t)  can  be  shown  to  be  the  unique  bounded  solution  to  the  initial 
value  heat  equation: 

u^  =  Au  ,  (1.3) 

u(x,0)  =  f(x)  .  (lo4) 

Next,  consider  v(x,t)  =  Au(x,t).   The   function  v  is   itself  a 
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solution   Co   the  heat  equation,  and  can  be  interpreted  In  three  wavs. 
First,  by  definition, 

v(x,t)  =  Au(x,t)  =  u^(x,t)  . 

Second,  by  standard  calculus  of  convolutions, 

v(x,t)  =  J    AK(x-y,t)f(y)  dy  .  (1.6) 

Here  we  think  of  v(x,t)  as  a  filtered  version  of  f(x),  using  Che  filter 
AK(x,t),  I.e.,  the  laplaclan  of  a  gausslan.  The  kernel  AK  can  be 
viewed  as  a  bandpass  filter, '■aAd'  Is  frequently  modeled  as  the 
difference  of  two  gausslans;  The  model  Is  logical,  since  K(x,t)  Is 
itself  a  solution  to  the  heat  equation,  and  thus 

AK(x,t)  =  K^(x,t)  5  l/e(K(x,t-K:)  -  K(x,t))  . 
Using  the  first  interpretation  of  v(x,t)  note  that 


-  J  v(x,t)dt  =  -  J  u^(x,t)dt  =  f(x)  (1.7) 

0  0 


provided  that  u(x,t)  ■>■    0  as  t  -►  » .   (It  is  sufficient,  for  example  that 
f  ^L^).  Thus  f(x)  can  be  reconstructed  from  v(x,t). 

The  second  Interpretation  views  v(x,t)  as  a  solution  to  the 
initial  value  heat  equation  using  initial  data  g(x)  =  Af(x).  We  see 
that  complete  knowledge  of  v(x,t)  is  assured  in  any  representation  from 
which  one  can  reconstruct  g(x).  Once  v(x,t)  is  determined,  then  f(x) 
(and  also  u(x,t))   are   also   specified,  according  to  Equation  (1.7). 
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Thus  a  representation  of  the  v(x,t)  data  provides  an  Implicit 
representation  of  the  f(x)  and  u(x,t)  data.  Since  v(x,t)  is  a  solution 
to  the  heat  equation,  we  are  henceforth  concerned  with  representations 
of  the  data  in  a  solution  to  the  heat  equation  such  that  some  form  of 
the  initial  data  can  be  reconstructed.  In  what  follows,  u(x,t)  will 
stand  for  any  generic  solution  to  the  heat  equation,  with  f(x)  as 
initial  data,  although  In  terms  of  the  above  notation,  v(x,t)  and  g(x) 
are  likely  instances  of  the  representation  problem.  ^ 

So  consider  a  solution  u(x,t)  to  the  heat  equation  defined  in  the 
space  {(x,t)  |xe  IR",  t  >  0} .  The  zero  set  of  u(x,t)  is  the  point  set  in 
(x,t)  where  u  =  0.  This  set  might  be  empty  (for  Instance,  if  f  Is 
subharmonic  or  superharmonic) ,  everything  (if  f(x)  Is  harmonic),  or  a 
proper  subset  of  (x,t)  space.  In  the  latter  case,  the  zeros  can  be 
Isolated  points,  lines,  and  surfaces  (but  never  regions).  However, 
components  of  the  zero  set  which  are  Isolated  points  and  lines  of 
codimenslon  greater  than  one  are  degenerate,  in  the  sense  that  u  will 
be  nonnegatlve  (or  nonpositive)  In  a  region  of  any  point  of  the 
component.  Components  of  the  zero  set  which  separate  positive  and 
negative  regions  will  be  surfaces,  and  have  a  special  definition: 

Definition: 

The  zero  crossings  refer  to  the  surfaces  separating  regions 
where  a  function  is  positive  and  where  it  Is  negative.   Foe  a 
solution  u(x,t)  to  the  heat  equation,  the  zero  crossings  refer  to 
the  point  set 

3{(x,t)|u(x,t)  >  0}  n  3{(x,t)|u(x,t)  <  0}. 
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Slnce  u(x,t)  is  a  smooth  function  for  t  >  0  (analytic,  in  fact),  the 
zero  crossings  will  necessarily  be  a  collection  of  surfaces  of 
codlinension  one,  analytic  except  possibly  for  discrete  cusps, 
satisfying  u  =  0  on  the  surfaces. 

For  imagery  data,  it  has  been  argued  that  the  zero  crossings  of  a 
laplacian  of  a  gausslan  filter  (AK(x,t))  of  the  image  data  Include  all 
edges  of  significance  at  the  corresponding  level  of  resolution  [4],  and 
accordingly  might  be  associated  with  human  visual  neural  processing  of 
retinal  intensity  data  [5].  Zero  crossings  have  also  been  used  as 
features  for  stereo  matching  and  motion  correspondence  between  pairs  of 
images  (e.g.,  see  [6]).  -j.-. 

It  is  clear  that  any  nonzero  scaj.ar  multiple  of  a  solution  u(x,t) 
will  have  an  identical  set  of  zero  crossings  as  u(x,t).  It  has  been 
recognized,  at  least  intuitively,  that  zero-crossings  alone  do  not 
contain  all  the  information  In  the  Initial  data,  even  modulo  a  scalar 
factor,  but  that  zero  crossings  together  with  slope  Information  should 
permit  reconstruction  [7].  Lack  of  a  theoretical  foundation  has  been 
lamented.  We  show  here  that  these  Intuitions  were  correct,  and  that 
the  theory  is  quite  simple. 

If  we  consider  a  representation  to  be  a  map  from  one  topological 
space  to  a  (hopefully  simpler)  topological  space,  the  value  of  the 
representation  can  be  determined  by  analyzing  the  fibers  of  the  map  and 
the  preimage  of  local  neighborhoods.  A  fiber  of  the  map  is  the  set  of 
elements  which  map  to  a  single  element  or  representation.  Generally,  a 
useful  representation  will  map  many  elements  to  a  single  element,  but 
that  the  class  of  elements  in  a  fiber  share  properties  which  are 
essential  to  the  interpretation  of  those  objects.   The  stability  of  the 
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representatlon,  however,  is  a  separate  matter.  Since  the 
representation  of  an  element  can  never  be  presented  with  infinite 
precision,  the  preiraage  of  a  local  neighborhood  of  a  representation  is 
the  true  collection  of  elements  represented.  Both  issues  need  to  be 
analyzed  for  the  proposal  to  represent  signal  and  image  data  by 
zero-crossings  of  bandpass  filters  of  the  data.  Our  results  in  the 
following  sections  are  intended  as  initial  investigations  along  these 
lines. 
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II«   Consequences  of  the  Maximum  Prlnctple 

In  this  section,  we  show  that  the  maximum  principle  for  parabolic 
equations  implies  that  zero  crossings  are  never  created  as  t  increases. 
We  further  show,  using  the  Hopf  maximum  principle,  that  there  exists  at 
most  one  solution  of  the  heat  equation  with  specified  zero  crossings 
and  specified  gradient  data  at  those  zero  crossings.  The  first  result 
has  been  discussed  in  the  context  of  data  representation  in  artificial 
intelligence  [8]  (at  least  for  n  =  1),  but  to  our  knowledge  previous 
treatments  were  independent  of  the  maximum  principle. 

Let  u(x,t),  xSlR^,  t  >  0  be  a  nonzero  bounded  solution  to  heat 
equation,  and  let  A  denote  its  zero  crossings.  The  zero  crossing 
surfaces  will  evolve  as  t  increases,  but  according  to  the  theorem 
below,  no  surface  can  be  completely  separated  from  the  Initial 
hyperplane  {t  =  0}: 

Theorem  2 . 1 :  A  component  of  A  is  a  surface  which  starts  at  the 
hyperplane  {t  =  0},  and  either  ends  at  the  hyperplane  {t  =  0}  or 
extends  to  infinity. 

Note  that  for  n  =  1,  the  zero  crossings  are  curves.  In  general, 
the  surfaces  of  A  admit  local  real  analytic  parameterlzations,  and 
therefore  have  at  most  isolated  cusp  singularities  .  The  content  of 
Theorem  2.1  can  be  summarized  as  3A  C  (t  =  0},  Before  proving  the 
theorem,  we  state  our  second  result: 

Theorem  2.2;  Suppose  we  are  given  a  collection  of  surfaces  A 
(consistent  with  Theorem  2.1)  and  specified  gradient  data  on  A.   Then 
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there   exists  at  most  one  solution  u(x,t)  to  the  heat  equation  having  A 
as  its  zero  crossings  and  the  specified  gradient  data  on  A. 

Thus  given  the  zero  crossings  and  gradient  data  on  those  zero 
crossings  of  a  solution  u(x,t)  to  the  heat  equation,  then  u(x,t)  is 
uniquely  specified,  and  there  is  no  other  solution  having 
simultaneously  the  same  zero  crossing  surfaces  and  gradient  data  on 
those  surfaces.  In  particular,  the  Initial  data  f(x)  =  u(x,0)  is 
uniquely  specified  by  the  zero  crossing  and  gradient  data.  The 
solution  u(x,t),  provided  it  exists,  can  be  constructed  by  the 
algorithm  presented  in  Section  III. 

Theorem  2.1  and  2.2  follow  easily  from  versions  of  the  maximum 
principle  for  parabolic  partial  differential  equations.  For 
completeness,  we  state  these  theorems  in  a  general  form  useful  for  us: 

Theorem  2.3  [as  stated  in  9]: 

Let  D  be  a  domain  of  the  (n+l)-dlmenslonal  (x,t)  space,  and 
suppose  that  Au  -  Uj^  >  0  in  D.  Suppose  that  the  maximum  of  u(x,t)  in  D 
is  obtained  at  an  interior  point  (xQ,tQ)  €  D,  with  uCxQ.tg)  =  M.  Then 
u(xpt^)  s  M  at  all  points  (xj.t^)  e  D  which  can  be  connected  to 
(xQjtg)  by  a  path  In  D  consisting  only  of  horizontal  and  vertical 
segments. 

Theorem  2.4.  A  Hopf  boundary  lemma  for  parabolic  equations  [9].  As  in 
Theorem  2.3,  suppose  that  Au  -  u^  >  0  in  a  domain  D.  Suppose  that  the 
maximum  M  of  u(x,t)  in  D  is  attained  at  a  point  P  €  8D.  Assume  that  an 
(n+1)  dimensional  sphere  In  (x,t)  space  can  be  constructed  such  that  P 
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lies  on  the  boundary  of  the  sphere,  and  Che  Interior  of  the  sphere  lies 

entirely   In  D,  with  u  <  M  In  the  sphere.   Also  suppose  that  the  radial 

direction  from  the  center  of  the  sphere  to  P  is   not   parallel   to   the 

t-axis.   Then  If  - —  denotes  any  directional  derivative  in  an  outward 

3v 

direction   on  3D,    we   have 


ll>  0     at     P. 
3v 


Proof  of  Theorem  2.1:  If  a  component  of  the  zero  crossing  surfaces  A 
encloses  a  region  D  separated  from  the  hyperplane  {t  =0},  then  D  is  an 
open  domain  with  u  =  0  on  3D.  Either  u  or  -u  (also  a  solution  to  the 
heat  equation)  will  contain  a  positive  maximum  attained  at  a  point 
(xg,tQ)  in  the  Interior  of  D.  Since  (xQ,tg)  can  be  connected  to  some 
point  of  a  bounding  zero  crossing  surface  by  a  path  whose  Interior  lies 
in  D,  consisting  of  vertical  and  horizontal  segments,  thus  Theorem  2.3 
is  violated.  • 

Proof  of  Theorem  2.2:  Suppose  u(x,t)  and  v(x,t)  are  distinct  solutions 
to  the  heat  equation  having  identical  zero  crossing  sets  A  and 
identical  gradients  on  A.  Then  w(x,t)  =  u(x, t)-v(x, t)  is  also  a 
solution  to  heat  equation,  vanishes  on  A,  and  has  zero  gradients  on  A. 
Now,  the  zeros  of  w  are  discrete  surfaces  and  points  which  include  A. 
Further,  A  consists  of  real  analytic  surfaces  with  possibly  discrete 
cusps.  Thus  there  exists  some  Pe  A  so  that  a  sphere  can  be  constructed 
through  P  lying  entirely  on  one  side  of  the  surface  of  A  on  which  P 
lies,  chosen  sufficiently  small  so  that  w  is  strictly  positive  (or 
strictly  negative)  inside  the  sphere,  and  finally  such  that  the   center 
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of   the  ball  is  not  directly  above  or  below  P.  Since  w  =  0  at  P,  either 

w  or  -w  has  a  maximum  In  the  closed  sphere  at  P.  By  Theorem  2.3,  _ — *    0 

a  V 

at  P,  in  any  outward  direction.   This  contradicts  the  fact  that  w(x,t) 
has  zero  gradient  data  at  P.  • 

All  of  the  theorems  in  this  section  hold  for  a  more  general   class 
of  parabolic  initial  value  problems,  i.e., 


Uj.  =  Lu  (2.1) 

u(x,0)  =  f(x)  (2.2) 

where 

"  32     "  3 

L  =  1     a..(x,t)  5 5 +1      b^(x,t)  5_ 

...   ^1      3x.3x4   .,   J      3x. 

is   a  strongly  elliptic  operator  with  bounded  coefficients.   In  analogy 
with  Equation  (1.1),  the  solution  to  (2. A),  (2.5)  can  be  written  as 


u(x,t)  =  i    G(x,t;y,0)f(y)dy  (2.3) 


where  G(x,t;y,T),  t  <  t  is  the  Green's  function  of  the  parabolic 
equation  (2.4).  The  heat  equation  (1.3)  is  special  in  that  its  Green's 
function  K(x-y,t-T)  (see  (1.2))  depends  only  on  I x-y |  and  t^ .  Thus 
the  solution  to  the  heat  equation  problem  (1.3),  (1.4)  is  given  by  a 
radially  symmetric  convolution,  in  agreement  with  our  intuitive  notions 
of  blurring  and  representation  at  lower  resolution. 
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III.   ReoonstruGtion  from  the  Zero  Croaslnga  and  Gradient  Data. 

Suppose  we  are  given  the  zero  crossing  set  A,  and  the  gradient 
data  for  a  solution  to  the  heat  equation.  In  this  section,  we  show  how 
to  reconstruct  the  solution  u(x,t). 

First,  consider  a  component  of  A  which  encloses  a  bounded  region  D 
of  the  (x,t)  half-space.  The  boundary  of  D  consists  of  zero  crossings 
r  C  A,  and  a  portion  Z  of  the  {t=0}  hyperplane:  3D  =  P  U  E.  Fix 
y  e  B '^  and  t  >  0,  with  t  greater  than  the  range  of  t  in  D.  Let 

v(x,t)  =  K(x-y,T-t),  (3.1) 

and  observe  that 


Vj.  +  Av  =  0   in  D.  (3.2) 


We  use  the  heat  equation  for  u  and  Green's  theorem  to  write 


0  =  J]  (u^-Au)v  +  (vj.-Hiv)u  dxdt 
D 

=  i   [  (uv)C  +  (uVv  -  v7u)«njda,  (3.3) 

3D 


where  do  is  surface  measure,  and  (n,C)(x,t)   is   the  outward  pointing 

g 
unit   normal  vector  to  3D  at  the  point  (x,t)  S  3D,  and  7  =  ,^  refers  to 

3  x 

spatial  derivatives  only.   Using  the  facts  that  (n  ,C  )  =  (0,-1)  on  Z   and 
u  =  0  on  r ,  equation  (3.3)  becomes 


i       u(x,0)v(x,0)dx  -  -  i  v(x,t)(VuTi)(x,t)da  (3.4) 

(x,0)€  z  r 

Let  us  define 


-13- 
f(x)    (x,o)e  z 

f(x)  =  {  (3.5) 

0       (x,0)^  Z 


and   u(x,t)  =  J    K(x-z, t )f (z)dz  .  (3.6) 

Then  equation  (3.4)  can  be  rewritten  as 


u(y,T)=J  K(x-y,  T-t)(VuTi)(x,t)da  (3.7) 

r 


Thus  if  we  know  T  and  the  gradient  of  u  along  T,  we  can  compute  u(x,t) 
for  t  >  max{t|(x,t)  G  D}.  Since  u(x,t)  and  f(x)  are  related  by  (3.6), 
we  are  in  a  situation  where  we  want  to  reconstruct  f(x)  from  u(x,t). 
This  is  achieved  by  solving  the  backwards  heat  equation.  John  [10], 
among  many  others,  have  considered  algorithms  to  solve  this  problem. 
In  general,  a  backsolution  does  not  exist  unless  certain  growth 
conditions  are  satisfied;  in  this  case,  the  existence  of  a  solution  Is 
guaranteed  since  f(x),  given  by  (3.5),  exists.  From  a  numerical  point 
of  view,  the  construction  of  f(x)  from  u(x,t)  is  unstable.  However, 
providing  f(x)  is  sufficiently  well  behaved,  since  f(x)  has  compact 
support  and  is  of  constant  sign,  reasonable  approximations  to  f(x)  are 
In  practice  possible. 

Once  f(x)  is  determined,  then  u  is  known  on  3D.  Theoretically, 
then,  u(x,t)  can  be  determined  in  D,  and  then  in  all  of  (x,t)-space  by 
analytic  continuation.  In  practice,  however,  closer  approximations 
will  be  obtained  by  treating  each  component  bounded  by  A  separately. 
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Since  Vu(x,t)  on  a  zero  crossing  surface  Is  always  normal  to  the 
surface,  the  gradient  data  can  be  specified  by  a  scalar  function  on  A. 
According  to  equation  (3.7),  the  most  convenient  form  for  the  gradient 
data  is  Vu-ndo,  where  7u  =  (-JL, . . .  ,_iL) ,  and  n  is  the  horizontal 
component  of  the  outward  normal  (n,C).  As  an  example,  consider  the 
special  case  n=l,  and  T  given  by  t  =  t(x),  |x|  <  a.  In  that  case, 
(3.7)  becomes 


a 

u(y,T)  =  j   K(x-y,  T-t(x))  iii(x,t(x))t'(x)dx  . 

3x 
-a 


Note  that  this  integral  can  be  written  in  the  form 


u(y,T)  =  ]  K(x-y,  T-t)  -li  (x,t)dt  , 
r  ^"^ 


where  the  Integral  is  over  all  (x,t)  ST.   A  similar  formula  holds   in 
higher  dimensions. 

However,  the  formula  (3.7)  applies  only  to  the  case  where  D  is 
bounded.  Next  consider  an  unbounded  component  D  of  (x,t)-space  with 
the  zero  crossing  surfaces  A  removed.  Let  t  be  greater  than  the  range 
of  t  in  all  of  the  bounded  components  of  (x,t)  space  less  A.  Finally, 
let  D^  =  D  n  {(x,t)|x  e  ]R",  0  <  t  <  t},  see  Figure  1. 


-15- 
We  once  again  consider  the  Integral  In   (3.3),   this   time   integrating 
over  D^ .   The   boundary  of   this   region  consists   of  zero  crossing 
surfaces  P,  portions  Z   C  {t=0},   and  an  additional   boundary   portion 
Z^  C  {t=T}.   Analogous  to  (3.4),  we  obtain 


i  u(x,0)v(x,0)dx  +  j  v(x,t)(VuTi)(x,t)d0 

(x,o)  ez  r 


u(y,T)     (y,T)  e  l^ 

(3.8) 
0         (y.T)  9'E^  . 


Here  we  have  used  the  fact  that  v(x,t)  =  K(x-y,0)  is  a  delta  function 
with  mass  at  x  =  y.  Define  f(x)  and  u(x,t)  as  in  (3.5),  where  Z  is  now 
the  {t=0}  portion  of  the  boundary  of  the  unbounded  component  D.  In  that 
case,  (3.8)  says 

u(y,T)  =j  K(y-x,  T-t)(7u«ri  )(x,  t)da ,   for  (y,T)$r  E^        (3.9) 

Since  u(y,T)  Is  unknown,  (3.8)  gives  no  useful  information  when 
(y,T)SE^^  To  obtain  the  u(y,T)  data  for  (y,T)e  E^  ,  we  must  content 
ourselves  with  analytic  continuation  of  the  u(y,T)  for  (y,T)f^E^, 
assisted  by  the  fact  that  equation  (3.9)  holds  with  the  unknown  but 
small  error  u(y,T)  when  (y,T)  e  E^ .  Once  u(y,T)  has  been  extended  in 
this  manner  for  y  €  IR",  reconstruction  of  f(x)  is  done  by  solving  the 
backwards  heat  problem. 

We  are  currently  conducting  numerical  experiments   to  test   the 
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feaslbility  and  accuracy  of  these  reconstructions  for  both  bounded  and 
unbounded  regions  enclosed  by  zero  crossing  surfaces. 
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IV.   Representation  using  Zero  Crossings  Alone. 

In  this  section,  we  sketch  some  Ideas  about  the  Inverse  problem 
that  arises  when  the  given  data  consists  of  the  locations  of  the  zero 
crossings,  with  no  Information  about  the  gradient  data.  Obviously,  the 
zero  crossings  of  a  function  represents  a  class  of  functions  that 
includes  all  scalar  multiples  of  the  original  function.  Thus 
reconstruction  In  this  section  refers  to  restoration  to  within  a 
multiplicative  constant. 

We  note  that  the  zero  sets  of  a  function  determine  the  function 
only  in  very  special  cases,  i.e.,  if  one  has  a  very  rich  analytic 
structure  such  as  complex  analyticlty.  For  example,  a  polynomial  is 
determined  (up  to  a  multiplicative  constant)  by  its  degree  and  full  set 
of  zeros.  Under  certain  conditions,  the  same  is  true  for  entire 
functions  in  the  complex  plane.  In  a  more  general  situation,  one  has 
the  celebrated  Riemann  Roch  theorem  of  algebraic  geometry.  Logan's 
theorem  [11)  also  is  a  uniqueness  theorem  in  a  certain  complex  analytic 
case.  We  emphasize  that  such  uniqueness  results  are  rare,  and  do  not 
hold  in  classes  of  solutions  to  general  elliptic  or  parabolic 
equations. 

As  a  simple  example  of  non-uniqueness  for  zero  crossings  of  the 
heat  equation,  consider  the  case  n  =»  I ,  A  =  {(x,t)|x  =  0,  t  >  0}. 
Clearly,  any  function  f(x)  which  is  odd  about  the  origin,  i.e., 
f(-x)  =  -f(x),  will  have  A  as  a  zero  crossing.  As  long  as  0  is  the 
unique  zero  of  f(x),  A  will  be  the  sole  zero  crossing  of  the 
corresponding  u(x,t).  It  is  not  hard  to  show,  using  Taylor  series 
expansion  of  f(x)  about  x  =•  0,  that  the  only  smooth  functions  yielding 
A  as  a  zero  crossing  are  odd  functions. 
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As   a   second  example  which  suggest  more  general  non-uniqueness, 
consider  a  solution  to  the  heat  equation,  u(x,t)  with  analytic   initial 
data   f(x),   n  =  1,   having  a   single  zero  crossing  curve  which  can  be 
parameterized  by 

t  =  t(x),  |x|  <  a  , 

where  t(a)  =  t(-a)  =  0,  and  t(x)  >  0  elsewhere.  We  further  suppose 
that  the  zero  crossing  curve  contains  no  cusps,  so  that  t(x)  is 
analytic.   We  have  that 

u(x,t(x))  =  0  for   |x|  <  a.  (4.1) 

Further,  since  u(x,t(x))  Is  an  analytic  function  of  x  for   |x|  <  a,   a 

necessary  and  sufficient  condition  for   (A.l)   to  hold  is  that  all 

d" 
derivatives  (u(x,t(x))   vanish  at  x  =  a. 

dx" 

For  simplicity,  we  assume   that   t'(a)  =  a   with  0  <  a  <  ».   The 

cases  a  =  0  and  a    =  <*>   can  be  treated  similarly.   Using  the  equations 

^t  "^  ^xx>  ^^'^  u(x,0)  =  f(x),  we  obtain  from  (4.1)  by  successive 
differentiations 


u^  +  u^*  t'(x)  =  0  , 
^xx+-^tx*t'  +  "tt*(t')2  =  0  , 


etc. 


Evaluating  at  x  =  a,  we  obtain  constraints  on  f(x); 
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f (a)  +  af"(a)  =  0 


f"(a)  +  of'"(a)  +  a2f""(a)  =  0 


Thus  f"(a)  is  determined  from  f (a),  and  f""(a)  is  determined  from 
f'"(a)  and  f (a).  In  general,  by  differentiating  n  times  and  evaluating 
at  X  =  a,  a  linear  constraint  of  the  form 

f(")(a)  +...+  o"f(2n)(3)  =  0  (4.2) 

will  be  obtained.  By  induction  we  see  that  the  even  derivative 
f(2n)(3)  j^g  determined  by  the  odd  derivatives  at  a  of  order  less  than 
2n. 

Now,  the  system  of  equations  (4.2)  is  also  sufficient  for  an 
analytic  function  f(x)  to  give  rise  to  a  solution  u(x,t)  satisfying 
(4.1).  Thus  any  analytic  function  satisfying  (4.2)  will  have  t  =  t(x) 
as  a  zero  crossing  curve.  Since  all  odd  derivatives  of  f(x)  at  x  =  a 
are  free  parameters  in  (4.2),  it  is  reasonable  to  expect  that  there 
exists  many  such  analytic  functions  f(x).  Of  course,  most  of  these 
solutions  will  have  other  zero  crossings.  However,  because  of  the 
large  number  of  free  parameters,  perhaps  perturbations  of  the  original 
expansion  for  f(x)  are  possible  which  yield  no  new  zero  crossings. 

We  are  currently  conducting  experiments  to  explore  degeneracies  in 
these  and  more  complicated  zero  crossing  surfaces  situations. 

Finally,  we  note  how  one  might  use  the  results  from  Section  III  to 
obtain  multiple  solutions  to  the  inverse  problem  obtained  by  specifying 
the  zero  crossing  surfaces  of  a  solution  to  the  heat  equation. 
Ideally,   one  would  like  to  augment   the  zero  crossing  data  with 
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arbltrary  gradient  data,  and  then  reconstruct  using  the  algorithm  of 
Section  III.  However,  this  procedure  will  In  general  fall  to  construct 
a  solution  with  the  same  zero  crossings  as  the  original  solution,  and 
may  even  fall  to  yield  data  which  can  be  backsolved  to  a  function  at 
t  =  0.  Instead,  the  gradient  data  must  be  constrained  to  satisfy  an 
Integral  equation. 

We  Illustrate  the  procedure  and  the  integral  condition  by  a  simple 
example.  Consider  the  case  n  =  1  with  a  single  smooth  zero  crossing 
surface  T  parameterized  by  x  =  x(t),  t  >  0,  starting  at  0,  x(0)  =  0  and 
extending  to".  The  curve  separates  the  (x,t)  domain  into  two  regions, 
Dj  on  the  left  and  D2  on  the  right  (see  Figure  2). 

We  will  determine  a  compatible  pair  of  data 

f(x)  =  u(x,0),  X  <  0  (4.3) 

Q(t)  =  u^(x(t),t),   t  >   0  (4.4) 

so  that  the  solution  to  the  heat  equation  defined  by 

t 
u(x,t)  =  J  K(x-x(s),  t-s)Q(s)ds 
0 

0 
+  i   K(x-y,t)f(y)dy  ,  (4.5) 

which  automatically  satisfies  (4.3),  will  also  vanish  on  T   and  satisfy 
(4.4). 

The  function  f(x),  x  <  0,  should  be  chosen  among  smooth  functions 
vanishing  at  0  and  -" ,  and  nonzero  elsewhere.  However,  f(x)  cannot  be 
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chosen  arbitrarily.   Instead,  it  must  be  chosen  so  that  there  exists   a 
solution  Q(t)  to  the  Volterra-type  Integral  equation 

t 
Q(t)/2  =  j  K^(x(t)  -  x(s),  t-s)Q(s)ds 
0 

0 
+  i   K(x(t)  -  y,t)f' (y)dy  (A. 6) 

In  that  case,  we  can  show  that  u(x,t)  defined  in  (4.5)  has  the  required 
zero  crossing  T,  and,  incidentally,  gradient  data  0(t).  Then  by 
choosing  different  compatible  pairs  f(x),  x  <  0,  and  Q(t),  different 
solutions  u(x,t)  can  be  obtained  each  giving  T  as  a  zero  crossing 
surface.  Of  course,  the  ability  to  find  one  such  pair  depends  on  the 
curve  r.  Not  all  curves  will  be  the  zero  crossing  belonging  to  some 
Initial  data. 
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Figure  1.   An  unbounded  component 
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